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$8 要 : 以 宁夏 石嘴山 市 引 黄 灌 区 为 研究 对 象 ,基于 灌 


区 内 76 眼 长 序列 地 下 水 监测 


,利用 经 典 统 计 与 地 统计 相 


结合 的 方法 ,对 灌区 地 下 水 埋 深 的 时 空 分 布 进行 分 析 , 并 通过 偏 相 关 分 析 确 定 其 影响 因素 。 结 果 表 明 :(1) 近 20a 
地 下 水 埋 深 整体 呈 增 大 趋势 , 增 大 总 幅度 为 0.09 m, 可 分 为 2 个 阶段 :1997 一 2002 年 为 逐渐 减 小 趋势 ,年 均 减 幅 为 
0.037 ma';2002 一 2017 年 为 增加 趋势 ,年 均 增幅 为 0.019 m-a; (2) 地 下 水 埋 深 空间 上 为 扇形 分 布 ,以 大 武 口 区 为 
原点 ,黄河 为 局 形 弧 段 ,从 西 往 东 逐渐 减 小 ,各 月 时 空 分 布 变 化 主要 受 区 域 灌溉 水 平 影 响 ; (3) 地 下 水 埋 深 空间 结构 
均 为 高 斯 模型 , 块 金 系 数 分 布 在 1.39%~3.08%, 为 强烈 的 空间 自 相关 性 ,空间 异 质 性 主要 受 灌区 本 身 的 地 质 结构 以 
及 补给 量 和 排泄 量 等 因素 影响 ,而 影响 程度 依次 为 :农业 灌溉 取 黄 河水 量 > 降雨 量 > 蒜 发 量 > 人 工 开采 。 
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宁夏 地 处 西北 和 干旱 缺 水 地 区 , 随 着 沿 黄 生 态 经 
济 带 的 快速 发 展 ,黄河 水 资源 日 趋 紧 张 ,干旱 缺 水 
问题 已 经 严重 制约 引 黄 灌区 经 济 和 社会 的 发 展 。 
石嘴山 市 是 银川 平原 降雨 量 最 少 .蒸发 最 强烈 Ek 
渍 化 程度 最 严重 地 区 ,目前 , 盐 渍 化 耕地 面积 4.06x 
10* hm , 占 总 耕地 面积 的 48.35%0m。 引 黄 灌 区 土壤 
盐 渍 化 根本 原因 是 区 工地 下 水 位 较 高 ”2 ,因此 , 控 
制 地 下 水 位 在 合理 的 深度 是 防治 土壤 盐 渍 化 发 生 
的 重要 途径 。 同 时 , 随 着 工 .农业 快速 发 展 .产业 
转型 升级 及 老 工业 基地 改造 加 快 ,水 资源 需求 量 剧 
增 , 地 下 水 资源 的 开发 利用 强度 逐步 扩大 "" ,宁夏 5 
个 地 下 水 超 采 区 中 石嘴山 分 布 4 个 钻 , 区 域 地 下 水 
超 采 问题 突出 。 

地 下 水 埋 深 受气 候 条 件 及 人 类 活动 影响 ,其 
中 ,降雨 量 "" 蒸发 量 灌溉 水 量 人工 开采 
被 认为 是 引 黄 灌区 地 下 水 埋 深 的 主要 影响 因素 ,但 
不 同 影响 因素 的 影响 程度 随 研究 区 域 不 同 而 不 
同 。 目 前 ,针对 石嘴山 市 引 黄 灌区 地 下 水 埋 深 的 影 
响 因 素 研 究 报道 较 少 ,因此 ,以 石嘴山 市 引 黄 灌区 
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76 眼 长 序列 地 下 水 监测 并 数据 为 基础 ,研究 地 下 
水 埋 深 时 空 演变 规律 及 其 影响 因素 ,以 期 为 该 区 
域 土壤 盐 渍 化 防治 .地 下 水 开发 利用 与 保护 提供 
utei m, 


材料 与 方法 


1.1 研究 区 概况 

石嘴山 市 位 于 宁夏 北部 ,地 处 105"58' ~ 
106°39'E,38°21' ~ 39"25'N ,为 我 国 西北 内 陆 区 , 东 
跨 黄 河 相 邻 内 蒙古 , 西 中 贺兰山 东 草 ,地 势 西高 东 
低 。 包 括 贺 兰 山区 平原 区 及 鄂尔多斯 高 原 区 ,其 中 ， 
平原 区 包括 山 前 洪 积 倾斜 平原 和 黄河 冲 湖 积 平原 。 
属 中 温带 干旱 气候 ,多 年 平均 降水 量 170.42 mm ,年 
内 降水 主要 集中 在 6—9 月 , 占 全 年 总 降水 量 的 
76.78%。 年 均 蒸发 量 为 1160.59 mm。 石 嘴 山 市 耕 
地 面积 8.38x10' hn? ,其 中 盐 渍 化 面积 4.06x10': hm’, 
占 比 48.35% ,石嘴山 市 2017 年 盐 涡 化 面积 统计 见 
表 1。 

根据 文献 [8] 可 知 ,石嘴山 市 2017 年 总 取水 
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表 1 2017 年 石嘴山 市 土壤 盐 渍 化 统计 
Tab.1 Statistics of soil salinization in Shizuishan City in 2017 


非 盐 化 盐 化 轻 盐 化 中 盐 化 重 盐 化 小 计 

面积 /hm 43287.54 40529.50 27363.08 9914.21 3252.21 83817.03 

比例 /% 51.65 48.35 32.65 11.83 3.88 100 
10.235x10* m’, 其 中 地 下 水 取水 1.047x10 mm , 占 总 上游 地 下 水 的 侧 向 补给 ,大气 降 水 及 引 黄河 水 灌溉 


取水 量 的 11.23%( 表 2)。 
1.2 水 文 地 质 概况 
石嘴山 市 位 于 银川 平原 断 陷 盆 地 ,第 四 系 分 布 
广泛 , 受 基底 构造 影响 ,第 四 系 沉 积 厚 度 呈 现 向 平 
原 边缘 变 注 的 规律 。 研 究 区 自 西向 东 可 分 为 山 前 
单一 潜水 区 和 中 东部 “ 双 层 结构 ”区 (图 1)。 单 一 洪 
apr a a JE ERINE URP E 
\ 砂 类 和 黏 性 土 组 成 ,含水 层 渗透 性 强 ,地 形 坡 度 
地 下 水 交替 强烈 ,主要 接受 基 岩 裂 际 水 .山区 沟 


BL 


直入 渗 补 给 等 。 
1.3 数据 来 源 

数据 由 宁夏 回族 自治 区 水 文 水 资源 勘测 局 提 
供 。 研 究 区 内 共 布 设 76 眼 长 序列 地 下 水 监测 井 (图 
2) ,每 个 月 用 皮 尺 与 测 绳 测定 6 次 (1 日 .5 日 、11 
日 .16 日 .21 日 .26 日), 每 次 重复 观测 3 次 ,每 月 地 


下 水 埋 深 值 为 6 次 观测 值 的 平均 值 。 时 间 分 布 及 影 
响 因 素 分 析 采 用 1997 一 2017 年 以 月 为 步 长 数据 , 空 
间 分 布 采用 2017 年 以 月 为 步 长 数据 。 


谷 常 年 性 地 表 水 流 及 洪流 的 垂直 渗入 补给 。 中 东 1.4 分 析 方 法 
部 "“ 双 层 结构 ?平原 区 地 势 平 坦 , 含 水 层 岩 性 颗粒 较 利用 SPSS 19.0 对 地 下 水 埋 深 进行 描述 性 统计 
细 , 渗 透 性 较 差 ,地 下 水 径流 缓慢 ,补给 来 源 主要 有 特征 及 地 下 水 埋 深 影响 因素 的 偏 相 关 分 析 。 半 方 
表 2 2017 年 石嘴山 市 取水 量 统计 
Tab.2 Statistics of water intake of Shizuishan in 2017 /105 n 
bh 农业 工业 城镇 生活 KHAA 生态 总 取水 量 
合计 “地 下 水 ”合计 地 下 水 合计 ”地 下 水 合计 ”地 下 水 (湖泊 补水 ) ”合计 ”地 下 水 
大 武 口 区 034 0.007 0.129 0.115 0.196 0.188 — 0.005 0.005 0.257 0.927 0315 
平 罗 县 6.577 0.0804 0209 0.164 0.07 0.05 0.035 0.035 0.04 6.931 0.333 
惠 农 区 1.763 0.000 0.527 0.213 0.072 0.072 0.015 0.015 2.377 0.399 
小 计 8.68 0.19 0.865 0.492 0.338 0.31 0.055 0.055 0.297 10.235 1.047 
地 下 水 取水 量 占 比 /% 2.19 56.88 91.72 100 10.23 
农 牧场 小 
<— aak ie 一 
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Fig. 1 Hydrogeological profile 
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A Fig.4 Variation characteristics of shallow groundwater 
depth in different months in Shizuishan City 


图 2 石嘴山 市 地 下 水 观测 井 分 布 示意 图 


Fig. 2 Distribution of groundwater observation well in 


Shizuishan City 


差 函 数 分 析 在 GS+9.0 中 完成 ,在 ArcGIS 10.0 中 Geo- 
Statistic Analyst 模 块 完成 克 里 格 插 值 。 


2 结果 与 分 析 


2.4 地 下 水 埋 深 时 间 变 化 特征 

2.1.1 年 际 变化 特征 ”根据 石嘴山 市 76 眼 地 下 水 
监测 井 长 序列 观测 数据 ,1997 一 2017 年 年 均 地 下 水 
埋 深 具有 明显 的 年 际 变化 特征 ,整体 为 增 大 趋势 ， 
增 大 总 幅度 为 0.09 m。 整 个 时 间 序 列 可 分 为 2 个 阶 
段 ,其 中 1997 一 2002 年 呈 逐 渐 减 小 趋势 , 减 小 幅度 
约 0.22 m ,年 均 变 幅 0.037 ma';2002 年 至 今 呈 逐渐 
增加 趋势 ,从 1.26 m(2002 年 ) 增 加 至 1.57 m(2017 
年 ), 总 增幅 0.31 m, 年 均 增幅 0.019 m-a (IL 3) 5 


年 份 
ASSI RNSINNENSINEEN AENEIS ESSRSERNYRNSNU RA RCRNARERNS RA 
YYRESRNNNS TENEI ANNANN 
二 AAAAAAANNANAAN 
1.10 y-0.01x--1.289 
g 120 R2—0.430 
È 1.30 
加 1.40 
1.50 
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图 3 1997 一 2017 年 石嘴山 市 浅 层 地 下 水 埋 深 年 际 变化 
Fig. 3 Annual variation characteristics of shallow groundwater 
depth in Shizuishan City from 1997 to 2017 


2.1.2 年 内 变化 特征 根据 各 监测 井 1997 一 2017 
年 各 月 地 下 水 埋 深 均值 变化 可 以 看 出 (图 4) ,在 整 


个 年 度 中 ,2 月 埋 深 最 大 ,为 2.05 mm, 之 后 随 气 温 逐 
渐 回 升 ,土壤 解冻 , 蔡 发 能 力 加 强 , 埋 深 逐 渐 减 小 ， 
但 变 幅 较 小 ,由 2.05 m 减 小 到 1.70 m; 从 4 月 下 旬 开 
始 春灌 , 埋 深 呈 急 速 减 小 趋势 ,5 月 达到 1.11 m, 较 4 
月 减 小 了 0.59 m, 之 后 6 一 8 月 达到 全 年 最 小 值 , 保 
TEE 1 m 左 右 ;9 月 田间 灌水 减少 埋 深 开始 增 大 ,10 
月 初 停止 放水 , 埋 深 增加 到 1.63 m; 从 10 月 底 开 始 
冬 灌 , 使 11 月 埋 深 再 1 次 迅速 减 小 ,达到 0.99 m, $7 
10 月 上 升 了 0.64 m. 
2.2 地 下 水 埋 深 描述 性 统计 特征 

由 表 3 可 知 , 各 年 埋 深 最 小 值 均 为 0, 最 大 值 在 
10.43 ~ 12.39 mo 2013—2017 年 埋 深 的 均值 分 别 为 
1.71 m,1.74 m,1.68 m,1.755 m 及 1.75m。 变 异 系数 
大 小 反映 总 体内 部 各 样本 之 间 的 变异 程度 ,一 般 分 
为 3 级 :Cv<10% 为 弱 变 异性 ,10%<Cv<100% 为 中 等 
变异 性 ,Cv>100% 为 强 变 异性 。 各 年 变异 系数 分 布 
在 3.89% ~4.96% , 均 为 弱 变 异性 ,说 明石 嘴 山 市 近 
年 地 下 水 埋 深 波动 较 小 。 

由 表 4 可 知 ,灌区 1 一 5 月 .9 月 和 10 月 变异 系数 
在 56% ~95%, 属 于 中 等 变异 性 ;6 一 8 月 和 11 一 12 
月 变异 系数 大 于 100% ,属于 强 变异 性 。 
2.3 地 下 水 埋 深 空间 结构 特征 

按照 决定 系数 较 大 、 残 差 较 小 的 原则 ,在 GS+ 
9.0 软 件 中 进行 半 方 差 吨 数 分 析 ( 表 5)。 块 金 值 Co 
为 随机 变异 量 , 块 金 系数 Cy(CotC) 表 示 由 随机 因素 
引起 的 空间 变异 占 总 变异 的 比例 , 即 空间 变异 程 
度 。 块 金 系 数 <25% 说 明 变 量具 有 强烈 的 空间 相关 
性 ,25% ~50% 说 明 变 量具 有 明显 的 空间 自 相 关 ， 
5096 ~75% 时 变量 具有 中 等 空间 自 相 关 , > 75% 时 
变量 空间 自 相关 性 微弱 ,变异 主要 由 随机 变异 组 
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表 3 石嘴山 市 地 下 水 埋 深 描述 性 统计 特征 


Tab. 3 Descriptive statistical characteristics of groundwater depth in Shizuishan City 


年 份 最 小 值 /mm ”最 大 值 /m 均值 /m 中 位 数 /m 方差 标准 差 变异 系数 /% 偏 度 
2013 年 0 10.43 1.71 1.53 1.29 76 3.89 20.26 
2014 年 0 11.51 1.74 1.5 1.45 84 4.26 2227 
2015 年 0 12.39 1.63 1.4 1.51 92 4.88 29.12 
2016 年 0 12.28 1.75 1.54 1.55 88 4.8 27.48 
2017 4F 0 12.72 1.75 1.52 1.54 88 4.96 29.68 


表 4 石嘴山 市 地 下 水 埋 深 各 月 描述 性 统计 特征 


Tab.4 Descriptive statistical characteristics of groundwater depth in each month in Shizuishan City 


月 份 最 小 值 /m ”最 大 值 /lm ”均值 /m 中 位 数 /m 方差 标准 差 变异 系数 /% 偏 度 I HE 
1H 0.7 12.18 1.73 1.8 1.68 1.29 65 5.26 32.39 
2H 0.66 12.24 2.05 2.1 1.62 1.27 56 5.19 32.12 
3H 0.47 12.28 1.95 203 1.7 1.3 60 4.99 30.72 
4H 0.33 12.38 1.7 1.76 1.92 1.39 70 4.81 28.12 
5H 0 12.36 1.11 1.23 2.14 1.46 95 4.86 27.51 
6 月 0 12.39 1.03 1.12 2.32 1.52 110 4.64 25.78 
7H 0 12.32 1.05 1.07 2.31 1.52 117 4.86 27.36 
8H 0 12.72 1.03 1.08 2.31 1:52 115 4.89 27.88 
9 H 0 12.68 1.38 1.52 2.11 1.45 83 4.81 27.93 
10H 0.58 12.7 1.63 1.86 1.92 1.38 66 4.87 28.8 
11 H 0.05 12.54 0.99 1.01 231 1.52 113 4.83 26.98 
12 H 0 12.33 1.26 1.21 2.16 1.47 101 5.08 29.44 
表 5 地 下 水 埋 深 半 方差 结构 
Tab.S Semi variance structure of groundwater depth 

月 份 模型 Co CotC 变 程 /km CoCotC/% R RSS 

| Gaussian 0.12 6.249 0.635 1.92 0.604 14.4 

2 Gaussian 0.12 6.249 0.654 1.92 0.595 13.8 

3 Gaussian 0.1 6.209 0.649 1.61 0.591 14.4 

4 Gaussian 0.11 7.229 0.692 1.52 0.605 15.1 

5 Gaussian 0.21 7.429 0.672 2.83 0.617 15.9 

6 Gaussian 0.19 7.389 0.644 2.57 0.628 174 

Y Gaussian 0.14 7.289 0.630 1.92 0.634 17.8 

8 Gaussian 0.2 7.409 0.640 2.70 0.616 18.4 

9 Gaussian 0.18 7.369 0.679 2.44 0.616 15.5 

10 Gaussian 0.15 6.309 0.653 2.38 0.607 13.4 

11 Gaussian 0:23 7.469 0.649 3.08 0.618 £7.7 

12 Gaussian 0.19 7.389 0.646 2.57 0.612 18.6 
年 均 Gaussian 0.1 7.209 0.670 1.39 0.62 15.6 


成 ,不 适合 采用 空间 插值 的 方法 进行 空间 预测 ”…。 ” 素 在 空间 上 的 变化 ) 引 起 。 

由 表 5 可 知 , 各 月 及 年 均 地 下 水 埋 深 均 符合 高 斯 模 2.4 地 下 水 埋 深 空间 变化 特征 

型 , 块 金 系数 分 布 在 1.39%~3.08%, 表 明 地 下 水 埋 深 ”2.4.1 年 均 空 间 分 布 特征 ”根据 确定 的 最 佳 变异 函 
具有 强烈 的 空间 相关 性 ,空间 异 质 主要 由 自 相关 因  ” 数 模 型 ,在 ArcGIS 软件 平台 通过 克 里 格 插值 ,对 
素 ( 灌 域 本 身 的 地 质 结构 以 及 补给 量 和 排泄 量 等 因 ”2017 年 石嘴山 市 年 均 地 下 水 埋 深 空间 分 布 特征 进 
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行 分 析 ( 图 5)。 由 图 5 可知, 地 下 水 埋 深 在 空间 上 整 
体 为 扇形 分 布 , 其 中 大 武 口 区 为 原点 ,黄河 为 扇形 
弧 段 ,由 西向 东 逐 渐 减 小 。>3 m 区 域 主要 分 布 在 
贺兰山 山 前 洪 积 书 ,大武 口 区 ,是 石嘴山 市 工业 聚 
集 区 ,也 是 地 下 水 漏斗 所 在 地 ;1~2 mm 区 域 主要 分 布 
在 平 罗 县 和 惠 农 区 引 黄 灌溉 区 , 约 为 总 面积 的 2/3， 
主要 受 黄河 水 灌溉 的 入 渗 补 给 ,地 下 水 长 期 处 于 高 
水 位 状态 ;2~3 m 区 域 分 布 在 >3 m. 1-2 m 区 域 之 
间 , 为 条 形 分 布 ;<1 m 区 域 主要 分 布 在 河 (黄河 ) , 湖 
( 沙 湖 等 ) 区 。 
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5 2017 年 年 均 地 下 水 埋 深 空间 分 布 
Fig. 5 The spatial distribution of annual average groundwater 


depth in 2017 


242 各 月 空间 分 布 特征 ”由 图 6 可 知 ,1 一 3 月 随 
温度 上 升 ,蒸发 能 力 增强 ,2~3 m 面积 逐渐 扩大 ;4 月 
随 农田 灌溉 期 的 到 来 , 埋 深 逐渐 减 小 ,1~2 m 面积 增 
大 ,2~3 m 面积 随 之 减 小 ;5 一 8 月 为 <1 m 和 1~2 mH 
深 面 积 交替 出 现 , 作 为 作物 主要 生长 期 , 随 灌溉 水 
量 的 累计 增加 ,<1m 埋 深 面 积 逐 渐 扩 大 ,7 月 达到 峰 
值 ,8 月 逐渐 降低 ,9 月 降低 到 与 5 月 相当 ,10 月 基本 
消失 ;从 10 月 20 日 冬 灌 开 始 ,11 月 埋 深 <l m 面积 
次 迅速 扩大 ;12 月 以 1~2 m 埋 深 面 积 为 主 , <1m 埋 
深 面 积 零星 分 布 ; 埋 深 >3 mm 区域 主要 集中 在 大 武 
口 区 附近 ,分 布 空间 相对 稳定 , 且 其 分 布 在 不 同月 
份 变 化 不 大 。 
2.5 地 下 水 埋 深 影响 因素 分 析 

石嘴山 市 引 黄 灌区 地 下 水 动态 变化 受降 水 28 
发 . 引 黄 灌溉 和 人 为 开采 等 因素 综合 影响 ,根据 地 


下 水 埋 深 与 降雨 量 KEE .农业 灌溉 取 黄 河水 量 、 
人 工 开采 量 的 偏 相关 分 析 ,确定 灌区 地 下 水 埋 深 影 
响 因 素 。 在 偏 相 关 分 析 中 , 目标 值 为 地 下 水 埋 深 ， 
自 变 量 分 别 为 降雨 量 .蒸发 量 .农业 灌溉 取 黄 河水 
量 及 人 工 开 采 量 ( 表 6) ,由 表 6 可 以 看 出 ,在 剔除 其 
他 变量 影响 条 件 下 ,确定 地 下 水 埋 深 与 降雨 量 L2 
发 量 .农业 灌溉 取 黄 河水 量 及 人 工 开 采 量 的 偏 相关 
系数 依次 为 -0.361 .0.228 ,-0.661 及 0.127, 即 影响 程 
度 依次 为 :农业 灌溉 取 黄 河水 量 > 降雨 量 > 蒸发 
E > 人 工 开采 ,其 中 农业 灌溉 取 黄 河水 量 与 降雨 量 
达到 0.01 显著 水 平 。 


3 ”讨论 与 结论 


地 下 水 是 区 域 重要 水 资源 ,也 是 影响 灌区 生态 
的 关键 因素 ,分 析 其 时 空 分 布 特征 及 影响 因素 ,对 
于 明确 区 域 地 下 水 可 持续 开发 利用 ,维持 生态 地 下 
水 位 具有 重要 意义 。 得 出 灌区 近年 地 下 水 埋 深 呈 
稳步 增 大 趋势 ,与 席 文 娟 ”在 此 区 域 研究 结果 一 
致 ,分 析 原 因 主 要 是 因为 近年 来 灌区 大 力 发 展 节 水 
灌溉 措施 导致 黄河 水 和 人 渗 补 给 减少 所 致 王 。 而 年 
内 地 下 水 埋 深 变化 及 各 月 地 下 水 埋 深 的 描述 性 统 
计 结 果 表 明年 内 变化 特征 符合 灌溉 型 ,此 规律 进 一 
步 证 明 灌 溉 水 量 是 影响 地 下 水 埋 深 的 第 一 要 素 。 
空间 分 布 结果 表明 ,灌区 地 下 水 埋 深 由 西向 东 逐 渐 
减 小 。 其 中 >3 m 地 下 水 埋 深 主要 分 布 在 贺兰山 山 
前 洪 积 扇 倾 斜 平原 , 大 武 口 区 周围 ,靠近 山 前 范围 
为 单一 潜水 结构 ,向 东 逐 渐 过 渡 到 潜水 OK HK 
层 结构 ,含水 层 富 水 性 强 , 水 质 良 好 ,易于 开采 ,所 
以 该 区 域 为 石嘴山 市 现 用 水 源 地 及 工业 用 地 下 水 
主要 集中 区 ,经 长 期 开采 已 形成 地 下 水 漏斗 ,同时 ， 
由 于 其 主要 受 基 岩 裂 际 水 .山区 沟谷 常年 性 地 表 水 
流 及 洪流 的 垂直 渗 人 补给 , 受 黄河 水 灌溉 入 活 补 给 
较 少 , 故 其 分 布 在 不 同月 份 变化 不 大 ,该 区 域 下 一 
步 应 尽快 完成 生活 、 工 业 用 地 下 水 替换 水 源 ,减少 
地 下 水 开采 ,以 解决 地 下 水 漏斗 问题 。1~2m 地 下 
水 埋 深 主要 分 布 在 黄河 冲 湖 积 平原 ,为 潜水 、 承 压 
水 多 层 结构 区 ,是 石嘴山 市 主要 农业 灌溉 区 , 受 大 
量 黄河 水 灌溉 入 渗 补 给 ,各 月 分 布 受 灌 溉 季节 影 
响 ,地 下 水 长 期 处 于 高 水 位 状态 ,土壤 盐 渍 化 问题 
突出 , 故 应 加 大 地 下 水 开采 力度 ,发 展 井 渠 结 合 江 
溉 ,减少 黄河 灌溉 水 量 。 本 文 就 灌区 地 下 水 资源 管 
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图 6 2017 年 1 一 12 月 月 均 地 下 水 埋 深 空间 分 布 
Fig. 6 The spatial distribution of average groundwater depth from January to December in 2017 
表 6 地 下 水 埋 深 影响 因素 偏 相 关 分 析 
Tab.6 The partial correlation analysis of influencing factors of groundwater depth 
控制 变量 影响 变量 相关 系数 显著 性 
RRE ,农业 治 溉 取 黄 河水 量 ` 人 工 开采 降雨 量 -0.361 0.005 
降雨 量 \ 农 业 灌 溉 取 黄 河水 量 、 人 工 开 采 蒸发 量 0.228 0.085 
降雨 量 、 蒸 发 量 、 人 工 开采 农业 灌溉 取 黄 河水 量 -0.661 0 
降雨 量 ,蒸发 量 .农业 灌溉 取 黄 河水 量 人 工 开 采 0.127 0.18 
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Temporal and spatial characteristics of groundwater and its influencing factors 
in the Yellow River diversion irrigation area of Ningxia Shizuishan 
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Abstract: Taking the Yellow River diversion irrigation area in Shizuishan City, Ningxia Province, as the research 
object, based on 76 long-term groundwater monitoring wells, this paper combines classic statistics and geostatistics 
to analyze the spatial- temporal distribution of the depth to groundwater in the irrigation area and determine its 
influencing factors by partial correlation analysis. The results indicate that (1) the depth of groundwater near 20 a 
shows an enlargement trend as a whole, with a total raise of 0.09 m, which can be divided into two stages: from 1997 
to 2002, it gradually decreased, with an average annual decrease of 0.037 m+ a`; from 2002 to 2017, it demonstrated 
an upward trend, with an average annual increase of 0. 019 m» a`'. (2) The groundwater depth is fan-shaped. With 
Dawukou District as the origin and the Yellow River as a fan-shaped arc, it gradually decreases from west to east. 
The spatial and temporal distribution of each month is mainly affected by the regional irrigation level. (3) The 
spatial structure of the depth to groundwater is a Gaussian model, and the nugget coefficient is distributed from 
1. 3996 to 3. 0896, which is a strong spatial autocorrelation. Moreover, spatial heterogeneity is mainly affected by 
the irrigation area' s geological structure andrecharge and excretion factors. The sequence of the impact degree is as 
follows: water taken from the Yellow River for agricultural irrigation? rainfall > evaporation > artificial exploitation. 
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